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Abstract 



o 

By integration of generalized BSBM and Brans-Dicke cosmological models, in this article, we 
investigate the theoretical framework of fine structure constant variation and current cosmic ac- 



celeration. We first develop a mathematical formalism to analyse the stability of the model. By 
employing observational data to constrain the model parameters, phase space is performed and the 
attractor solutions of the model are detected. We then examine the model against observational 
data such as observational Hubble parameter dataset and quasar absorption spectra. The results 
confirms current universe acceleration and also predicts fine structure constant variation. Further- 
more, extrapolation of the best fitted model in high redshift (z > 15) illustrates a significantly 
larger variation of fine structure constant in earlier epoch of the universe. 
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1. INTRODUCTION 



There have been many researches of the cosmological consequences of allowing some of 
the physical constants of nature varying in spacetime. One of the first ones is that of varying 
the Newtonian gravitation constant, G, which leads to the Brans-Dicke (BD) scalar-tensor 
theory of gravity A second quantity the electron charge, and/or the fine structure 



constant, a = e 2 /he, in cosmological time where originally introduced in {2]. In general, 
we have no enough information of why the constants take their current values, or whether 
they are truly constant or not. Nevertheless, in case of fine structure constant variation, 
motivated by recent accurate observational evidence from the quasar absorption spectra Si- 



ll], cosmologists seek to find models to justify the variation. Due to the first observational 



evidence, the fine structure constant might change with time; smaller than its present value 
by ^ = ~ 10~ 5 |l2| -|14l| . A varying a due to variation of electron charge originally 



proposed by Bekenstein 15j and then after the first observational evidence from quasar 
absorption spectra revived and generalized by Sandvik, Barrow and Magueijo in the so 
called BSBM proposal 16]. In BSBM model, an scalar field coupled to the electromagnetic 
part of the Lagrangian is responsible for the alteration of fine structure constant. Later, a 
generalized BSBM model has been studied by Barrow in jl?!] where the coupling constant u 
in the theory is assumed to be a function of the scalar field in the model. 

From a different approach, in 18], we investigated both universe acceleration and fine 
structure constant variation in BSBM theory in the presence of an exponential self potential 
for the scalar field in BSBM model. We performed stability analysis and constrained the 
model with the observational data. The attractor solutions are detected for the best fitted 
model parameters. We then checked the model against Hubble parameter dataset and quasar 
absorption spectra. The numerical computation verified both redshift dependence in the fine 
structure constant and also current cosmin acceleration. 

Barrow et al. [19] also in a different work investigated a cosmological model that in- 
corporates variation of both Newtonian gravitation constant and fine structure constant. 
Following [l9|, in a new framework in this manuscript, we also integrate both Brans-Dicke 
(BD) and generalized BSBM theories in which the gravitational constant G is replaced by 
the BD scalar field <fi and fine structure constant is related to BSBM scalar field i/j. In 



comparison with the work in 19], here, we assume a generalized coupling constant as a 



2 



function of scalar field u(i/j). In addition, we constrain the model parameters by using re 



cent observational data of Type la Supernovae 



20|- 22J and perform stability and phase 



space investigation. The stability analysis urges us to find attractor solutions of the dynam 



ical system as asymptotically stable solutions 



231 ]- 30] . We therefore avoid fine tuning issue 



in which the dynamical trajectories in the phase space lie along a common track despite 
starting from different initial conditions. We finally check the model against observational 



data of hubble parameter 



31| and quasar absorption spectra. 



2. STABILITY OF THE MODEL AND CONSTRAINTS ON THE PARAMETERS 

We start with the integration of Brans-Dicke and BSBM models. The action describing 
the dynamics of the system with varying both fine structure constant, a and gravitational 
constant, G, takes the form 

S = J d'xy^^R - ujbd^P + (1) 

where 

C = -^W^ + £ em e(- 2 ^ + C m . (2) 

The fine structure constant, a, and gravitational constant, G, respectively are given by two 
different scalar fields; a = aoe 2 ^ and G = 1/0. We assume that the parameter u(ip) = |f, 
as a function of scalar field ip, still determines the strength of coupling between scalar field 
and photons. The characteristic length scale, I, is introduced for dimensional reasons and 
gives the scale down to which the electric field around a point charge accurately obeys 
Coulomb force law. From the present experimental constraints, in order to avoid conflict 
with observation, the corresponding energy scale, y, has to lie between a few tens of MeV 
and the Planck scale ~ 10 19 GeV. A varying u(ip) is due to variation of length scale, I. In 
the conventional BSBM theory, variation in a is assigned to the electron change e. The 
negative sign of the first term in the Lagrangian density, (j2J), usually represents a phantom 
field. However, as we see later, our model does not produce a phantom regime in the history 
of the universe as there also exist other factors contributing to the dynamics of the system. 
The electromagnetic lagrangian is C em = —\f fJ , u f' J ' l/ where f^ u is the electromagnetic field 



3 



tensor. From action (JT]), the field equations in flat FRW cosmologies are 



3 F 2 v 1 6 2 ^ 



2H + 3H 2 = -X-- Pr e~^ - ^-V) ~ 2H^ - - (4) 



where we put 8tt = c = h = 1 and assume a perfect fluid filled the cosmic with p m = jp m . 
The ratio £ = ^y 2 ^ in equation ([2]) describes the fraction of nonrelativistic matter in the 
universe that contributes to C em . The energy density p m stands for the contribution from 
cold dark matter (CDM) to the energy density. We also neglect radiation energy density, 
p r = 0. In addition, the field equations for the scalar fields and ip respectively are 

» + 3 gj=»»< 1 + ™E . (5) 

3 + 2u BD 6 + 4UBD 

■ (du/dibU 2 ICI 9 ,„ 
2u u 

To reduce the complexity of the nonlinear field equations ([3])-(j6])and also perform stability 
analysis and plot phase space trajectories, we introduce the following new variables: 

v 2 = Pm P 2 = wj?_ 2 = n2 = t ( 7 \ 

X 30# 2 ' ? 0# 2 ' 1 2 # 2 ' u { ' 

Employing the Friedmann constraint equation fl3]), in terms of new variables, the fine struc- 
ture constant takes the from: 

a = x 2 «o|CI(l - X 2 ~ j + rj - cosn^y 1 . (8) 
For further consideration, in term of new variables the J^- is found to be, 

TH = -«X 1 + « + + ubdV + ^— ^ f + 277. 9 

In terms of these variables, the field equations dS])-© become, 

X = -§(»? + 3 + 2^), (10) 

C = -kv + S + 2§- 2 )-6^ X W, (11) 
2 H z a 

( , o , ^ , 3x 2 a + ^)-e 2 n9 , 

' = -^ + 3 + ^)+ 3 + 2^ ' (12) 

tf = o{n-p), (13) 
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where ' denotes derivative with respect to N = ln(a) and (3 is a dimensionless constant in 
oj oc e^ N . Next, we perform stability analysis while constraining the model parameters by 
using the distance modulus of 557 Type la supernovae (SNela) from the Union 2 sample. 
Note that ^ and in the above equations are given by relations flUJ) and (jSJ). Therefore, 
the only parameters in the model are ojbd, P, in addition to the initial conditions for 
the dynamical variables, \i £j 1i 6 an d Hubble parameter h . To derive constraints on 
cosmological parameters we also need the luminosity distance (g^) given as a function of 
the above model parameters and initial conditions. From numerical calculation, Table I 
shows the best-fitted model parameters, initial conditions and \ 2 minimum value. 



parameters 


^BD 


P 


x(o) 


£(o) 


?7(0) 


0(0) 


h 


v 2 • 

Aram 




-1.6 


-0.4 


0.4 


0.2 


-1 


±2 x 10" 6 


0.702 


543.8467 



Table I: Constraints on mode 



parameters from SNela data. 



The system, modeled by the autonomous differential equations ( ITDl - (I13j) . depends on 
the numerical values of the parameters that appear in the equations, we already constrained 
these parameters. By performing linear stability analysis, we find seven critical points given 
in Table II: 



Critical points 


(x, £, v, 0) 


PI 


(0, 0, 0, 0) 


P2 A ,b 


(0, ±3E, -6, 0) 


P3a,b 


(0, ±E, -2, 0) 


P±A,B 


(± F 1 0) 



Table II: Characteristics of the critical points 



^t-ubd + 3w£ D ). From Tables I and II, all the 



where E = ^/-2(3 + 2co BD ) and F = ^2 
critical points for the best fitted model parameters are real. The eigenvlaues for the first five 
critical points are respectively, (0, — (3, —3/2, —3/2), (6, —3, 3/2, —f3 — 6) and (1, 2, 3/2, —2 — 
0) where for the best fitted parameter /3 = —0.4, these points are unstable. On the other 
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hand, the eigenvalues for critical points P^a,b are (-^, — f^j^f, ^ 
stability conditions for these two symmetric critical points are (3 < 




where, for the best fitted parameters (3 and ujbd, are stable points. In a multidimensional 
phase space, a single critical point characterizes the entire system at an instant of time. 
The set of trajectories in phase space for all possible initial conditions and for a given set 
of stability parameters forms phase portrait of the system. From stability point of view, 
for attractors the trajectories rapidly return to the stable points after finite perturbations. 
The basin of attraction of an attractor is the set of all points for which the trajectories 
asymptotically approach it. For a physical interpretation, it is more appropriate to analyze 
the lower dimensional phase space which is formed by the intersection of phase portrait 
with an hypersurface. In our model, the 4-dim phase space is therefore reduced to a 3-dim 
one by projecting the phase space (x^ 7 ?^) into the subspace (77, £, %) at 9 = 0. In Fig.l, 
qualitatively the behavior of these solutions are exhibited by sketching the phase space 
portrait for the stable critical point. 



Fig 1: The projection of 4-dim phase space into 3-dim phase plane for the stable 
symmetric critical points P^a,b with the best fitted model parameter. 

The phase space trajectories usually define a set of unstable and stable (attractor) critical 
points. Following a small perturbation, the trajectories quickly return to the symmetric 
(with respect to x = line) attractors P^a,b- The red trajectory is the best fitted one 
with the observational data, leaving the unstable critical point P3a and entering the stable 
critical point PA a- A more detailed discussion of the physical properties of the system with 
cosmological application is given in the next section. 



best filed trajectory 
unstable point 
stable point 




6 



3. DISCUSSION 



In this section, for the model under consideration, we test the dynamical behavior of the 
cosmological quantities. The reconstructed effective equation of state parameter, w in terms 
of the phase space dynamical variables given by 



w 



i 



-2(1 + oW 



CKh 2 2 2 

2u B d ol 3 3(3 + 2w B£) j 



(14) 



From numerical calculation, the parameter w is shown in Fig.2. The best fitted trajectory 
(dashed red line) shows that the universe leaves the unstable state in very high redshift ( 
radiation dominated era) and enters the current quintessence stable state. The model does 
not exhibit phantom crossing behavior in the past, but predict that the universe eventually 
reaches the stable quintessence era at the present era (with wq = —0.66 within the limit of 



our observations estimation 



321]) . The other trajectories correspond to the dynamics of the 



universe for arbitrary (not fitted) model parameters. 




redshift 



Fig. 2: The reconstructed equation of state parameter, w , as a function of redshift. 
The best fitted parameter leaves the unstable state in the past corresponding to 
the critical point P%a and enters the current sable state corresponding to the 

stable critical point P^a 

In direct contact with observation the time shift density parameter — is the next quantity 
to be checked against observation. In terms of new variables, the parameter is given by, 



Act 

a 



ICIx 2 



[q 



'h 2 + 2r? 



ubdV z 



(15) 



2(3+2u BD ) 2' 

From numerical calculation, in Fig 3, the parameter is compared with its quasar absorption 



spectra observational dataset We see that the best fitted model parameter, while pass- 



ing through the datapoints, is not very sensitive to a variation in 10~ 5 scale. The model 
extrapolation predicts that for redshift z > 15 the time shift density parameter begins to 
vary significantly. Again, for comparison, we also plotted the fine structure constant for 
arbitrary (not fitted) model parameters. 




redshift 



Fig. 3: The time shift density parameter, Aa/a, plotted for the model against 
quasar absorption spectra observations. The best fitted trajectory 
is shown with red dashed line. 



Furthermore, we compare the behavior of Hubble 
calculation in the model with the observational data 
Hubble parameter in the model agrees with the data. 



parameter derived from numerical 



331 ] . From Fig. 4 the best fitted 



#160 




Fig. 4: The graph of derived Hubble parameter if as a function of redshift 
in comparison with the data. The best fitted parameters is shown with a red 

dashed curve. 

Finally, the reconstruct scalar field responsible for both a variation and universe acceler- 
ation using the best fitted model parameters is shown in Fig. 5. 
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W 1 10° 10 1 IO 2 ID 3 10 4 

Fig. 5: The graph of the scalar field ip as a function of redshift. The best fitted scalar field 
with the model parameters is shown with a red dashed curve. 

4. SUMMARY AND REMARKS 

In this paper we have introduced a criterion that can be used to probe the cosmological 
viability of gravitational and fine structure constants ( G, and a) variation theories. Specif- 
ically, we investigated the integration of BD and BSBM theories where the two BD and 
BSBM scalar fields are responsible both for the universe acceleration and a variation. In 
practice, the model has simultaneously satisfied both observational evidence from SNela and 
quasar absorption spectra. We first best fitted the model with the observational data from 
SNela for distance modulus and found constraints on the model parameters. We then per- 
formed two quantitative and qualitative analyses. The validity of model with the constraints 
on its parameters are verified by these analyses. Qualitatively, we performed stability anal- 
ysis and found the attractor solutions in the phase space. These attractor solutions start 
with a small perturbation of the model parameters in the radiation dominated epoch and 
exhibit the final resting state of the universe in an accelerated expansion phase. The matter 
dominated phase of the universe is just a transient state in this model. Two quantitative 
tests are also performed to directly compare the model with the observational data; the 
observational Hubble parameter and a variation tests. The best fitted model with the at- 
tractor property are verified by these two tests. The model also predicts that variation of 
fine structure constant become greater in higher redshift. From equation of state parameter 
the universe is in quintessence regime at the present time. This is in consistent with the 
attractor property of the equilibrium solution in a scalar field dominated cosmology. This is 
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an advantage over those cosmological models that predict universe acceleration in phantom 
era where the null energy is violated. In addition, by extrapolating backward in time the 
model predicts a significantly larger variation of fine structure constant in earlier epoch of 
the universe. 
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